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Abstract 

A new generalization of the t-J model with a nearest-neigbour 
hopping is formulated and solved exactly by the Bethe-ansatz method 
in the thermodymanic limit. The model describes the dymanics of 
fermions with different spins and with isotropic and anisotropic inter- 
actions. 
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In recent years there has been considerable interest in studying low- di- 
mensional electronic models of strong correlation due to the possilility that 
the normal state of the two-dimensional novel superconductivity may share 
some interesting features of a ID interacting electron system [1]. In one 
dimension, the Bethe-ansatz technique can allow one to exactly solve Hamil- 
tonians in special cases, such as the Hubbard model [2] and the ordinary 
t — J model at its supersymmetric point [3,4]. For example, it is possible to 
obtain the low-energy gapless excitation spectrum around the ground state 
by the finite-size scaling method [5,6] and calculate the critical exponents of 
the correlation functions [7-9]. 

The t — J model is a lattice model on the restricted electronic Hilbert 
space, where the occurrence of two electrons on the same lattice site is for- 
bidden. This restriction corresponds to an implicitly infinite on-site Coulomb 
repulsion. Two types of interactions between electrons on nearest-neighbours 
sites are considered: a charge interaction of strength V and a spin- exchange 
interaction J. The Hamiltonian of the extended version of the t — J model 
has the form [3,4,10] 
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where Cj )OL annihilates an electron with a spin component a, P is the projector 
on the subspace of non-doubly occupied states, and L is the lattice size. We 
have introduced an anisotropy in the charge interactions through a matrix 

In the isotropic case V a p = V the Hamiltonian (1) corresponds to the 
traditional t — J model which was exactly solved by the Bethe-ansatz method 
at the supersymmetric point (V = —J = 1) for the case S = 1/2. [3,4,9-11]. 
The generalization of this result for the arbitrary spin S was carried out in 
[12-14]. Other generalizations of the t — J model were studied in [15 -18]. In 
particular, in [18] the anisotropic generalization of the t — J model has been 
constructed and it was shown that the model (1) is solvable for the arbitrary 
spin and special values of the coupling J and V a p : 

J = e o, 
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V a p = e [(l + £ a) cosh 7 • <5 Q/3 + exp[s^n(a - /3)7](1 - <J a/? )] (2) 

where 7 > is a measure of the anisotropy. It was shown more exactly in [18] 
that the Hamiltonian (1,2) is the quantum counterpart of the so called Perk- 
Schultz model [19] which was diagonalized by Schultz in the most general 
form [20] (see also [21-22]). In fact, in paper [19] Perk and Schultz have 
considered two models, the Hamiltonian of the first one is given by (1,2) and 
the Hamiltonian of the second one is given by [19,22] 

L N 

H =~J2 2 4a c j{> c j+l,a c j+l,P> ( 3 ) 
j=l a,f3=l 

The systems which are described by (3) and their different modifications were 
studied in spin formulation in [23-25] (see also references therein). In this 
letter we present a new set of models of strongly-correlated particles which 
are exactly solvable. In the system there are different types of particles and 
interactions between two particles in the neighbour sites are given by either 
Hamiltonian (1) or Hamiltonian (3) and depend on the type of particles. 
Thus we will consider the system the Hamiltonian of which contains the 
interactions of both types (1) and (3) . 

Firstly let us formulate the problem and write down the Hamiltonian. 
Consider a periodic one-dimensional lattice of L sites. Place n particles 
on this lattice and specify its nature. We assume P components, denoted 
Qi, Q2, Qp and define riQ i as the number of particles of component Qi 

p 

(J2 n Qj = n )- 

i=i 

In each set there is Nq 1 

(J2N Ql = N). 
i=i 

types of particles, so that the summation 

n Qi +...+n q . 

means 

aeQi a=N Ql +...+N Q ._ 1 +l 
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We establish the one-to-one correspondence between particles a a(a,a G 
Qj). So in each set there are qQ j conjugate pairs of particles 

[(N Qi + l)/2] < q Qi < N Qi 

where [N/2] means the integer part of the number N/2. The dynamics of 
above defined particles is described by the Hamiltonian 

L N 

H = - E E P (cja^+l,« + 4+l,* C J>) p 
j=l a =l 

"EE ^ e i[ U a U p )c j a c jP c j+l,S C j+l^ 

~ E E e o{toc+ a c^c+ +1 ^c j+ i iQ -exp[s^n(a-/3)7]n ja n i+1)/3 }, 

(4) 

where £j, e = ±land the parameters C/W pl a y the role of anisotropies inside 
of each set and satisfy 

£tf> = l/E/£\ E(^) 2 = 2cosh 7 (5) 

The ordinary t — J model given in (1,2) is obtained by choosing sector of 
(4) in which there are particles only of one type from each set Qj and for this 
type a ^ a. The anisotropic version of the Perk-Schultz model of the second 
type (3) is obtained by choosing a sector in which we have particles only of 
one component, for example, Q\ and N = L. Previously we solved exactly 
the model (4) at this case also for N < L [25,26] for the different choice 
of conjugate pairs. The model which is described by the Hamiltonian (4) is 
also the generalization of the model considered by Sutherland [10] for the case 
when particles of the a given component are not identical and the anisotropy 
has been introduced in the system. In the general case the Hamiltonian (4) 
describes the dinamics of fermions with different spins. For example, the case 
Nq. = 2 corresponds to S — 1/2 , the case Nq. = 3 corresponds to S — 1 
with biquadratic interactions [26]. For the general case Si = (Nq — l)/2 the 
magnetic interactions inside of each set can be written as a polynomial of 
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degree 2Si in the spin operator. Moreover the interactions between particles 
of different sets also exist. 

The exact solution for the eigenstates and eigenvalues of the Hamiltonian 
(4) can be obtained within the framework of the Bethe-ansatz method [27,28] 
The central object of this method is the two-particle scattering matrix S 
which is calculated from the single- and two-particle processes described by 
the Hamiltonian (4). The nonvanishing elements of the S-matrix are 

S^ih, fc 2 ) = [sin(i7 - Ai + X 2 )]~ 1 S^,(X 1 - A 2 ); 

- SiS^Sp^'U^U^ sin A; for a, (3 e Qi 

Sa?/3'W = *<W'<W sinh 7 exp [isign((3 - a)X\ 

- G ik 5 aa >5pp> sin A; for a e Q { \ (3 e Q k ,i ^ k (6) 

where Xj (j = l,2,...,n) are suitable particle rapidities related to the mo- 
menta {kj} of the electrons by 

_ fn-G(Xj^i), e = -1, m 
^■-\-e(A i; j7), eo = +l, [n 

with the function G defined by 

0(A;7) = 2arctan(cot7 ■ tan A) ; -n < 6(A, 7) < it. (8) 

A necessary and sufficient condition for the applicability of the Bethe-ansatz 
method is the Yang-Baxter equations [27,29 ]. Up to our knowledge the 
form of the S-matrix (6) is a new one, therefore it is necessary to check 
the Yang-Baxter equations. We have checked these equations numerically 
for the different choice of parameters U$ and the sets Qi. We may note 
also that S-matrix (6) can be obtained by the baxterization of the following 
Hamiltonian 



JV-l 



H = E 



3 1 

ej = cosh 7 + £ ]T [G lk EfE^ 1 -smh 1 stgn(a-(3)E« a E^ 1 ], 

i^k a£Qi,f3£Q k 

+ E E efp^ufEfEf^-^E-E^- (9) 
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where the N x N matrices E a @ have elements 

(E al3 ) y s = 8 ai 5f3s 

The quantum chain (9) have important property that if ej satisfies the fol- 
lowing Hecke algebra [31] 

[e j: e k ] = 0, forjj - fc| > 2 

e) = 2 cosh 7^ (10) 

then the S-matrix (6) satisfies the Yang-Baxter equations and the model (4) 
is integrable [31]. We have checked that ej satisfies the equations (10) and 
thus we have proved the integrability of the model under consideration (4). 

The Hamiltonian (4) is diagonalized by a standard procedure by imposing 
periodic boundary conditions on the Bethe function. These boundary con- 
ditions can be expressed in terms of the transfer matrix of the non-uniform 
models (6) by using the quantum method of the inverse problem [31,32]. 

The rapidities Xj that define a n-particle wave function are obtained by 
solving the equations 



sinh(Aj — if/ 2) 



L 



(-l) n_1 A(Aj), (11) 



sinh(Aj + ry/2) 
where A (A) is the eigenvalue of the transfer matrix 

n 

r# } } (A) = E II 5 S +1 ( A < - A )> (/W = (12) 

It is simple to verify that besides the numbers of particles in each set n Qi 
the numbers of "conjugate" pairs in each set are conserved quantities in the 
Hamiltonian (4). Here we denote two conjugate particles of same set paired 
if they are consecutive particles and have no unpaired particles of this set 
between them. 

In a general case, the complete diagonalization of the transfer matrix (12) 
is not a simple problem even in the simplest special cases. Therefore here we 
restrict ourself to consideration of the thermodynamic limit of the Hamilto- 
nian (4), when the periodic and free boundary conditions are identical. 
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Firstly consider model (4) in the sector where we have no pairs of par- 
ticles. If we have no pairs of a set Qi then the first interaction term in the 
Hamiltonian (4) does not work and all particles of this set are identical and 
can be considered as one component of the model. In this way the general 
model (4) in this sector can be reduced to the anisotropic t — J model with 
P components [18,20] and the diagonalization of the transfer matrix of the 
inhomogeneous model (12) gives the following Bethe-ansatz equations 



x 



sm(A 3 - X-, + ^7) _ 

11 . 777?) TT^T) ■ r - ~ e o e ^+i 111WpW<t+iJ 

j>=i sm(X) - A}, - ^7) p=o 

- sinjxf } - + ie g+l7 ) "ft 1 sin(Af } - A^ +1) - f e CT+l7 ) . 

11 . /,((T) ,(<t) • , 11 . /,(<t) \(<t+1) . i \' ^ ^ 
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where 



j = l,2,...,m ff ; a = 0, 1. 
A ^ = 0; A^ = Xj] 

rij = rrij-i — rrij] m\—L\ m q = 0; ; n Q = n (14) 



and 



q — P; e — 1; G ffff — GW — G^o — 1; 

ei=£u G ik = g ik ; (i, k = 1, 2, P) (15) 

rij is the number of particles from the set Qj. 

The total energy and momentum of the model are given in terms of the 
particle rapidities Xj in the following form 

n n / sinh^ 'V \ 

E = —2 cos fc 7 - = 2g£i I cosh 7 — 



( , j=1 v cosh 7 - cos 2Aj / ' 



n 



P = ^(Xj). (16) 

Consider now the model (4) in the general case when in each sector Qj 
we have rij particles and n'j "conjugate" pairs of particles. The reference 
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state \l/o is made up of a state in which there are no "conjugate" pairs. 
Examining (4 ) we see that when this Hamiltonian acts on a state, it looks 
for "conjugate" pairs and replaces them by a sum over all such pairs from a 
given set Qj. The second possible process is the permutation of two neighbour 
particles or one particle and one conjugate pair. It is important that both 
processes do not depend on the number of types of particles in each set Nq 3 
up to the boundary. Thus the Bethe-ansatz equations depend only on rij 
and n'j .Therefore we may consider the Hamiltonian (4) for the case when 
all Nq j = 2. Certainly the fact that the solution in the general case is 
strongly in line with the solution for the case Nqj = 2 is connected with 
the Hecke algebra. The two models with open boundary conditions have 
the same energy levels with different multiplications if they have the same 
underlying Hecke algebra [30]. Thus we consider the Hamiltonian (4) in the 
thermodynamic limit when the periodic and open boundary conditions are 
equivalent for the case when all Nq = 2. The Bethe-anzatz equations for 
this case will be valid for the arbitrary choice of Nq . . The examinination of 
the case Nq. = 2 shows that the Hamiltonian ( 4) can be reduced directly 
to the 2P-component t — J model [18,20]. Thus the general solution of the 
model ( 4) is given by (13 -14) and (16) with redefined parameters e a and 

q = 2P; 6q — 1; Gcrcr = Go a = G a o — 1; 
ti = £[(i+i)/2]; G ik = #[(i+i)/2],[(fc+i)/2]; [(i + l)/2] ^ [(k + l)/2]; 

(i,fc = l,2,..,2P); 
G 2 i-i )2 j = G 2l ,2i-i = e,; (/ = 1,2,...,P) (17) 

and now nn-\ and n-u —n<n-\ are numbers of "conjugate" pairs and separate 
particles of the set Qi(l = 1, 2, P). 

The ground state energy and the excited states of the Hamiltonian (4) 
can be calculated in principle by straightforward methods on the base of 
Bethe-ansatz equations (13-14) and (17), which in the thermodynamic limit 
can be written as integral equations. However, the solution of these equtions 
depends strongly on sign functions e a and the phase factors and can be 
subject of separate considerations for the different choice of these parameters. 
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